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Abstract 

The structural and electronic properties in common-anion 
(GaSb)i/(InSb)i and common-cation (InAs)i/(InSb)i [111] ordered super- 
lattices have been determined using the local density total energy full 
potential linearized augmented plane wave method. The influence of 
the ordering direction, strain conditions and atomic substitution on the 
electronic properties of technological and experimental interest (such as 
energy band-gaps and charge carrier localization in the different sublattices) 
were determined. The results show an appreciable energy band-gap narrowing 
compared to the band-gap averaged over the constituent binaries, either in 
[001] ordered structures or (more markedly) in the [111] systems; moreover 
energy band-gaps show an increasing trend as the substrate lattice parameter 
is decreased. Finally, the systems examined offer interesting opportunities 
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for band-gap tuning as a function of the growth condition (about 0.7 eV in 
(GaSb)i/(InSb)i and 0.3 eV in (InAs)i/(InSb)i). 



I. INTRODUCTION 



Ternary systems based on III-V semiconductors (such as disordered alloys |T] [|. 



heterostructures |§-f| or quantum well systems ||10Hl2]1) have been the subject of wide 
scientific interest and of accurate theoretical studies, since they could be used as fundamental 
components in a large class of important devices (laser diodes or infrared detectors, to 
name just a few) |fL3| . In the present work, we focus our attention on superlattices 
(SL), whose structural, electronic and transport properties can be opportunely tuned 
by varying the constituent materials, the strain, the ordering direction or the layers 
thickness. To this end, we have examined the properties of interest in ultrathin [111] 
ordered SL, specifically in common-anion (GaSb)i/(InSb)! systems and in common- 
cation (InAs)i/(InSb)! systems, using the self-consistent all-electron FLAPW method 
within the density functional formalism. The systems considered here are now under 
experimental investigation (results obtained for strained-layer InSb/GaSb quantum well 
fLUl , In^Gai-^Sb/GaSb heterostructures fl5|] and InAsi.^Sba, alloys [|l6i|T7| have already 



been published); at the present time, however, we are not yet able to compare our predicted 
results with experimental values regarding the SL. 

The appreciable mismatch between the lattice parameters of the binary constituents 
(5.7 % in (GaSb)i/(InSb)i and 6.4 % in (InAs)i/(InSb)i) gives the opportunity to study 
the effects of the strain conditions on the SL electronic properties. In analogy with the 
common experimental approach [|T^ , p~8| , p~9[] we have considered various growth conditions for 



the SL, leading to different strain modes in the structure: (i) pseudomorphic growth on a 
substrate usually constituted by one of the binary constituents (in which the lattice constant 
parallel to the growth plane is taken equal to that of the bulk semiconductor composing the 
substrate) and (ii) a "free standing mode" (FSM), in which no constraints are imposed on the 
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bond-lengths, leading to relaxed lattice constants for the binary constituents, both different 
from their bulk values. The structural parameters for all of the structures considered have 
been chosen through total energy minimization or according to the macroscopic theory of 
elasticity (MTE) 0. 

The choice made towards the [111] direction is encouraged by recent experimental 
observations ^TJ that suggested the spontaneous ordering along this particular axis (the so 



called CuPt structure) shown by some III-V alloys during vapour-phase growth. The strong 
influence of the ordering direction on the SL electronic properties is immediately clear if we 
consider the ternary SL Brillouin Zones (BZ) as obtained from the binary zincblende BZ 
through folding operations, which are obviously different for [111] or [001] ordered systems; 
the immediate consequence of this is a noticeable difference in the electronic properties. 

In order to study the dependence of the electronic properties on the ordering direction, 
we have studied the common-anion and the common-cation systems in the three different 
growth conditions both in the CuAu-like (having [001] direction as growth axis) and CuPt- 
like structure. 

Following the model proposed by Wei and Zunger P%| , we consider the SL as obtained 



first from an ideal virtual crystal (a common-anion (common-cation) system having the 
cation (anion) with intermediate properties between the two cations relative to the binary 
constituents). We then introduce a perturbative potential, having a structural part (due to 
atomic displacements and relaxation due to epitaxial strain) and a chemical one (due to the 
electronegativity difference between the constituent atoms). 

In order to separate the effects due to the two different terms in the expression of the 
potential, we have studied strained (AC)i/ (AC)i-type systems ideally obtained by monolayer 
deposition of the same binary constituent (AC) along the [111] direction, in which the two 
different AC bond-lengths are equal to those in the equivalent (AC)i/(BC)i SL. 

We will discuss the results obtained in this work as follows: first of all, we will briefly 
expose the computational details and parameters used in the calculations (Section II); in 
Section III our structural results will be reported for all the different systems considered. In 



Section IV we will discuss the electronic properties of the SL, with particular attention to 
the quantities of technological interest (such as energy band-gaps and crystal field splittings) 
and their dependence on growth direction, atomic substitution and strain conditions. We 
will also discuss the charge density distribution and in particular the localization of the 
charge carriers in the different constituents sublattices. Section V summarizes our main 
results and draws some conclusions. 



II. METHOD OF CALCULATION 

We have determined the properties of the structures considered using the density 
functional formalism, within the local density approximation (LDA) |23| exchange and 
correlation potential as parametrized by Hedin-Lundqvist | 24}| . The calculations were 



performed using the ab initio all-electron full-potential linearized augmented plane wave 
method (FLAPW) |2^| . Core electrons as well as valence ones are treated using a self- 
consistent procedure; the shallow Ga 3d and In Ad states are considered as valence states, 
for which scalar-relativistic effects are included in the self-consistent calculation, whereas 
spin-orbit effects are treated in a perturbative approach. For [111] ordered systems, angular 
momenta up to l max = 6 in the muffin tin spheres (with radius R a = 2.4 a.u. for all the 
constituents atoms) and plane waves with wave vector up to k max = 3.3 a.u. are used, 
leading to about 600 basis functions. 

To perform integrations in reciprocal space, a set of four special k points is chosen in 
the trigonal Brillouin zone (BZ), following the Monkhorst-Pack scheme [26|| . Similar values 
for these computational parameters have been used for [001] ordered systems, with the only 
exceptions represented by l max = 8 and a set of three special k points used for the integration 
over the tetragonal Brillouin zone. Finally, the Broyden |?7[ method is used to accelerate 
the convergence in the self-consistent iterations. 
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III. STRUCTURAL PROPERTIES 



The atomic ordering along the [111] direction of a SL grown on a (111) substrate gives 
origin to a trigonal Bravais lattice with C$ v (Schoenflies notation) space group P8 |. The 



unit cell in real space contains 4 atoms and the origin is taken on a cation site [29 . 

The free structural parameters are determined following the macroscopic theory of 
elasticity (MTE) [f2T]], taking into account the elastic properties of the constituent materials, 
and then compared with those obtained through total energy minimization. We observe 
that in each cell, there are two atoms belonging to the same chemical species (the two 
Sb anions in the (GaSb)i/(InSb)i systems and the two In cations in the (InAs)i/(InSb)i 
systems), which are not equivalent from the coordination point of view. As an example, 
we consider the particular case of (GaSb) 1 /(InSb) 1 SL, having the first Sbc a bound with 
three Ga and one In, and the second Sb/ n showing a complementary situation. The total 
energy minimization procedure that considers all the free parameters in the unit cell as 
different degrees of freedom is a very onerous computational problem; this encouraged some 
simplifications, such as considering equal bond lengths between equal atomic species (i.e., 
in the common-anion system, we have chosen di n sb Ga = di n sb In ), thus reducing to three the 
number of degrees of freedom (in-plane lattice parameter and two different bond-lengths). 
Although this simplification is frequently used in total energy minimization PSj , one should 



be aware that this approximation results in considering an average over the two different 
local environments and therefore that the real elastic structures may slightly differ from our 
optimized SLs. 



According to MTE, using the same notation as in Ref. p0| , the structural parameters 
for the epilayer are determined as follows: 
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where Cjj are the elastic constants for the bulk epilayer (we have used the experimental 
values reported in Ref . |3(J ) . 

Through total energy minimization of the "ideal" (AC)i/(BC)i unrelaxed structures (in 
which all the atoms are arranged in a zincblende structure with lattice constant a - our free 
parameter - and with bond-lengths cLac — dsc — a (\/3/4) ), we have found an in-plane 
lattice constant very close to the average of the bulk constituents, according to Vegard's 
rule. We have thus examined a free standing mode structure (indicated in the following as 
Elastically Relaxed or simply ER), that has this value for the in-plane lattice constant. In 
Table [I] and Table [TI| we report the calculated structural parameters for the ternary common- 
anion and common-cation systems: the SI (S2) system is a common-anion SL grown on a 
GaSb (InSb) substrate, while the S3 (S4) system is a common-cation SL grown on an InAs 
(InSb) substrate. 

In the case of pseudomorphic growth on a substrate, we have found general agreement 
between the structures obtained through total energy minimization and those given by MTE; 
this fully justifies our having considered this approximation to determine the five unknown 
parameters in the free standing mode (FSM) structure. Note that deviations from the 
results predicted by MTE occur in the case of InSb strained to GaSb or to InAs; in both 
these structures, total energy minimization gives an InSb bond-length that is systematically 
larger (within 0.6 %) than the one expected according to MTE, even though the difference 
between the total energies for the elastic SL and for the total energy minimal structures is 
very small (barely larger than the numerical uncertainty of 1 mRy/unit cell). However, this 
can be justified considering that, due to its elastic properties, this material could easily be 
out of the linear elastic region. In fact, the elastic constants for InSb are quite smaller than 



those for GaSb and InAs [30], resulting in a larger effective strain - due to the mismatch - 



in the former case; this is also confirmed by the non-linear behaviour of the band-gap as a 



function of the strain, as will be discussed later. 



As expected from the similarity of the GaSb and InAs elastic constants [30] and bulk 



moduli (the experimental values are B InAs = 0.579 Mb and B GaSb = 0.578 Mb pl|), we 



obtain similar deformations for these two constituents respectively in the common-anion 
and common-cation SLs. We also notice that strains (either parallel or perpendicular) 
and percentage deviations from bulk bond-lengths are more pronounced in the common- 
cation systems, compared to the common-anion systems: this is obviously a consequence 
of the greater mismatch between the constituent lattice parameters in the (InAs)i/(InSb)i 
structures. 

In the case of [001] ordered systems (whose structural parameters are reported in Table | 
and Table |I| respectively for common-anion and common-cation SL), we obtain a tetragonal 
Bravais lattice with D\ d space group |28j and a unit cell in real space with 4 atoms (two 



of which are equivalent); the origin is taken on a cation (anion) site for the common-anion 
(common-cation) system. 

The MTE relations reported in Eqs.([I|) are still valid for [001] ordered SL, if the parameter 
D is redefined as: 

£,[ooi] = 2 / £i2\ (2) 
Vcn / 

A comparison between the [111] and [001] ordered structures (having the same chemical 
composition) shows that in the same growth conditions (FSM or growth on a substrate) the 
structural parameters are not equal. In particular we notice larger perpendicular strains in 
the [001] compared to the [111] ordered structures, while the parallel strains are obviously 
equal in considering the same growth conditions; what we find is thus a smaller deviation 
from bulk bond-lengths, due to a more effective relaxation. 



IV. ELECTRONIC PROPERTIES 



7 



A. Electronic levels 



The determination of the SL electronic energy levels is a fundamental point for most of 
the properties of interest in the systems considered. In Table [TTI| we report the calculated 
electronic levels (with a numerical uncertainty of ± 0.04 eV, equal for all the energies 
reported in the present work, unless otherwise specified) at the BZ center (T), for the 
different [111] ordered systems considered (free standing mode and pseudomorphic growth 
on the two substrates), both for the common-anion and for the common-cation systems. We 
also report the zincblende state from which the SL state derives through folding the f.c.c. 
Brillouin zone back into the smaller ternary trigonal zone. 

The splitting, Aqf, of the triply degenerate Ti 5v zincblende state is due to the non- 

(2) 

cubic crystal field and is conventionally taken positive if the doubly degenerate state T 3v 



has a higher energy compared to the state Ti v . As can be seen from Table we obtain 
a negative Aqf m the case of GaSb (InAs) strained to InSb for common-anion (cation) 
systems - corresponding to an in-plane extensive strain en - whereas the complementary case 
(pseudomorphic growth on a GaSb-substrate (InAs-substrate)) and the free standing mode 
produce a positive Aqf- 

(2) 

The introduction of spin-orbit coupling removes the double degeneracy of the T 3 J state 
and yields the electronic energy levels illustrated in Fig.|I] as a function of the substrate 
lattice constant. The topmost valence bands (E 1; E 2 , E 3 ) have been labeled according to 
the "quasi-cubic" model |E| (taking the centre of gravity of the SL valence bands as zero); 

Ei = +|(A S . D . + Acf) 

E 2i3 = -|(A S . . + Acf) ± §{(A S . . + A CF f - § A s . . A C f} 1/2 
considering Aqf as obtained from Table |TTT| and A s-0 . for the SL as the value averaged over 



the equivalent calculated p3| quantities for the binary constituents (Even neglecting the 



A s . . negative bowing, occurring in the common-cation SL ||22|| , it is possible to uniquely 
identify the SL levels with the "quasi-cubic" ones, which differ at most by 0.05 eV). From 
Fig.|l| we notice that the Ei and E c (the lowest conduction state) levels show an almost linear 

8 



behaviour as a function of the substrate lattice constant, while the other valence band states 
show a more complex trend, due to the interplay between crystal-field and s.o. effects (see 
the cross-over between the Ei and E2 states). 

In Table [TV] we report the relevant electronic energy levels for the CuAu systems at T; 
the notation is analogous to that of Table |TTT| , where folding relations of the f.c.c. Brillouin 
zone in the now tetragonal ternary zone involve the zincblende states as indicated in the 
Table. The trend in the signs of the crystal-field splittings Acf is similar to that evidenced 
in the [111] ordered systems (see Table [PV|) . 

The general underestimate of the band-gap energy in LDA has stimulated many attempts 
to solve this problem, but correction algorithms p4|-|36[| need an extraordinary computational 



effort in the SL case; thus our LDA band-gap energy (E^ A ) was corrected starting from the 



experimental data of the binary constituents | 37| . Due to a lack of experimental band-gaps 
regarding strained binaries, we have fitted the calculated values obtained for each binary 
in different strain conditions (i.e. tetragonal and trigonal), assuming a linear trend for 
the band-gap energy as a function of the in-plane strain ey: E^ A (e\\) = E^ A (0) + a ey, 
where Eg^ A (0) is the binary equilibrium calculated band-gap (a parabolic trend Eg^ A (e\\) = 
Egap A (®) + a e \\ + P e y has been used for InSb, which is assumed to be out of the linear 
region). Once we determined the coefficient a, we translated the curve so that it becomes: 
-^gap P ( e ||) = -^gap*(0) + a e lh w h ere -Egap*(0) i s ^e binary equilibrium experimental band- 
gap. We have thus used these empirical values to obtain the empirical band-gap energy 
averaged over the strained binaries (< E e ™^ >). Summing this quantity to the correction 
(5 = Eg^ A (SL)— < Eg® p A > ), we have finally obtained the predicted band-gap energy in 
the SL ( Eg™P(SL) ) with a numerical uncertainty of ± 0.05 eV. Although this procedure is 
empirical, it is expected to give good estimates of the real band-gaps, since it is well known 
that while the band-gap is strongly underestimated, the band-gap behaviour as a function 
of pressure is always very well reproduced by LDA JEJ . 

We report in Tables [V] and [VI| (respectively for [111] and [001] ordered structures), the 
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band-gap energies as obtained from LDA self-consistent unperturbed calculations (Eg£g), 
with the introduction of the perturbation due to spin-orbit coupling (Eg^ A ) and with the 
correction starting from experimental data (-E|™ p ). From these Tables we first notice that 
for all the systems considered we find a negative E g ^ A (due to an inversion which causes 
the conduction band minimum (CBM) to lie below the valence band maximum (VBM)). 
Furthermore, we observe that the larger the substrate lattice parameter, the smaller the 
band-gap energy (either in common-anion or in common-cation systems). 

A comparison between Table |V| and Table [Vl] confirms the trend predicted by Wei and 



Zunger []37fl for the energy band-gap: 

e%£ ] < < E Z ( 3 ) 

where E g 1 ^ and E g ^ are the band-gap energies respectively in the [111] and in the 
[001] ordered structures, while Eg™ is the band-gap average energy taken over the binary 
constituents (the calculated values of the LDA band-gap energy for the binary constituents 
are E^ A (GaSb) = - 0.47 eV, E^ A (InSb) = - 0.67 eV and E^ p A (InAs) = - 0.63 eV). 

As is well-known [^2J, band folding in the superstructures causes a repulsion between 
two binary electronic states of different symmetries, folded on a state of the same symmetry 
in the ternary phase and coupled through the perturbative potential mentioned in Section 
I (in Tables |T| and |V| the superscripts (1) and (2) indicate the two states involved in the 
repulsion mechanism). One of its interesting effects is the band-gap narrowing, compared to 
the equivalent quantity averaged over the binary constituents (as confirmed by the second 
inequality in Eq.([jD). The amount of this effect |22| is inversely proportional to the difference 



[e(ri c ) — e(Li c )] (in the [111] structure) or to the difference [e(ri c ) — e(Xi c )] (in the [001] 
structure); this difference is smaller in the [111] structure, causing a more striking band-gap 
narrowing than in the [001] structure (as shown by the first inequality in Eq.(|3|)). These 
observations are confirmed by the calculated values for the band-gap bowing parameters 
(defined, in analogy with the 50%-50% alloys, as b gap = A(E™p — E gap )) reported in Tables 
[V|and[V]]: we obtain a larger bowing in the [111] structures compared to the [001] ones and, 
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looking at the constituent chemical species, we can say that the bowing in common-cation 
systems is larger than in the common-anion ones. 

The band-gap trend as a function of the substrate lattice constant and its dependence 
on the ordering direction have been illustrated in Fig.|2], where we report the LDA band- 
gap [E LDA - solid line) and corrected band-gap (E emp - dashed line) as a function of 
the substrate lattice constant for (GaSb)i/(InSb)i (Fig.|| (a)) and (InAs)i/(InSb)i (Fig.|| 
(b)). A comparison between the SL energy band-gaps and the average value (E ave ) of the 
experimental (LDA) band-gap in the pure binaries - indicated by the filled (empty) circles - 
clearly shows the band- gap narrowing effect. 

Tables [V|and [VI] also show that the crystal field splittings in the [111] ordered structures 



are always bigger than in the [001] structures in the same growth conditions, with the only 
exception represented by the S4 system, which has a smaller Aqf compared to the other 
structures. This apparently strange behaviour can be explained by considering that the 
rfj state interacts with the lower T^j state, resulting in an upward shift (an effect relevant 
only in the common-cation SLs); furthermore, in the S4 system, the VBM is a r lt) state, 
which is only slightly involved in the level repulsion mechanism. Thus the stronger the level 

(2) 

repulsion, the larger the T^ v upward shift and the smaller Aqf becomes: this observation 
is thus a further proof of the validity of the band repulsion model. 

We observe that in the [111] ordered SLs the band-gap energy is determined by the 
difference in energy between the VBM - slightly localized on the anion belonging to the InSb 
monolayer (as will be clearly shown in the next section) - and the CBM - strongly localized 
on the GaSb (InAs) sublattice in the common-anion (cation) superlattice. Thus we could 
think of (GaSb)i/(GaSb)i-type systems as common-anion SLs in which we substitute the 
InSb monolayer with a GaSb monolayer. We would expect in this case a small modification 
of the band-gap energy, since the VBM will no longer be localized on the Sb belonging to 
the InSb sublattice but rather on the Sb belonging to a GaSb sublattice; therefore, this 
will be only a second order effect. Our prediction is confirmed by the calculated band-gap 
energies (spin-orbit included) reported in Table [VII| (second column) which prove the almost 
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total independence of the VBM on the cationic substitution; the change of the band-gap 
energy in the different structures is thus caused by the structural term in the perturbative 
SL potential rather than by the chemical term. An equivalent interpretation considers the 
(InSb)i/(InSb)i-type systems as SLs in which we have substituted the GaSb monolayer 
(where the CBM is strongly localized) with an InSb monolayer. In this case, the cationic 
substitution implies the chemical alteration of one of the atomic species (Ga) on which the 
wave function is strongly localized; thus, what we expect, is an appreciable change in the 
band-gap energy, as confirmed by the third column in Table |VII1] (from which we notice the 



E^R. increase). 

An analogous trend is observed for the common-cation systems, where InAs has now 
susbstituted the GaSb as the InSb partner in the SL (see second and third column in Table 

fm> . 

The trend in the crystal field splitting is strongly dependent on the class of systems 
considered. In fact, in the common-anion systems, this quantity is almost independent of 
the cationic substitution (as expected, because of the anionic character of the VBM, localized 
on the Sb atom). In this case, the chemical term of the potential existing in the SL has very 
little effect on A C f, compared to the structural term (see Table |VTJ). On the other hand, in 
common-cation systems anionic substitution has a strong effect on the crystal field splitting 
and the chemical term in the SL potential is now much more important than before, even 
though the structural term still has a strong effect on A C f (see Table |VIII|) . 



B. Charge density distribution 

One of the main effects of the perturbative potential in the SL (in particular of its 
chemical term, due to differences in the constituent atom's orbital energies ||) is the 
localization of the charge density in one of the constituent monolayers, which varies from 
state to state. As an obvious consequence, this effect causes the confinement of the charge 
carriers (holes or electrons respectively for valence or conduction states) in a different 
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sublattice. 

In order to better clarify the character relative to the different states of interest, their 
angular decomposition - for the common-anion systems in the three growth conditions 
considered - is reported in Table |D£| (we do not to report the equivalent Table for 
(InAs)i/(InSb)i, since this system is very similar to the previous one, as far as the charge 
decomposition is concerned). Referring to the charge density of the T\ v state, we notice in 
particular the growing s character on the InSb monolayer and the decreasing p character on 
the In atom as the substrate lattice parameter is increased; at the same time, the s charge 
density on the GaSb sublattice decreases, while the p charge grows on the Ga atom. 

We report in Fig.|3| (a) distribution for the Ti v state for the (GaSb)i/(InSb)i elastically 
relaxed (ER) systems, drawn the same for all the charge densities reported in this work in a 
plane perpendicular to the atomic layers. This state comes from p z orbitals (as we can see 
from the typical "butterfly" shape along the vertical growth ^-direction) and shows a strong 
bonding character, between different monolayers and within each monolayer. 

In FigJ^ (b) we report the charge density distribution relative to the T 3v (VBM) state for 
the common-anion system in its elastically relaxed structure, where the localization of the 
charge density in the InSb sublattice is particularly evident. We notice that this peculiarity is 
much more enhanced in the common-cation system (not shown), as a probable consequence 
of the anionic character of this state: in fact, what we expect in the common-cation system 
is for the Sb-atom to draw more charge than the As-atom. We have found that the charge 
density distribution in this state is not strongly influenced by the strain conditions, as can 



be seen from Table IX 



The calculated charge density distribution for the first conduction state r^, (CBM), 
relative to the common-anion elastically relaxed system, is presented in Fig.|] (a). What is 
relevant in this state is the strong localization of the charge density in the GaSb monolayer 
(a similar behaviour is shown by the common-cation ER system, where the charge density 
is concentrated on the InAs monolayer). 

The localization emphasized above becomes more pronounced as the substrate lattice 
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parameter is increased: the charge density distribution concentrates more and more in the 
GaSb monolayer (InAs monolayer) while at the same time the InSb monolayer becomes 
charge-depleted (as confirmed by Table [TXp . The second conduction state, r^r, shows 
a complementary trend, owing to the charge density distribution that is more and more 
concentrated on the InSb sublattice as the substrate lattice parameter is increased (as we 
notice from Fig.f| (b) for the common-anion elastically relaxed system; this behaviour is 
similar to the common-cation systems). 

As a consequence of these observations, in all these structures we have a direct gap in 
reciprocal space, while we obtain a "spatially indirect" gap, due to the localization of the 
rjjj state (VBM) on the InSb sublattice and of the r£? state (CBM) on the GaSb (InAs) 
sublattice in the common-anion (common-cation) SL. 



V. CONCLUSIONS 

Ab initio FLAPW calculations, based on density functional theory within LDA, have 
been performed in order to determine the electronic properties of ultrathin SLs. In 
particular we have studied a common-anion (GaSb) \j (InSb) i system and a common-cation 
(InAs) 1/ (InSb) i system, ordered along two different ([111] and [001]) directions . 

The relevant results obtained for these structures can be summarized as follows: 

1. Both the (GaSb)i/(InSb)i and (InAs) 1/ (InSb)! systems show a direct gap (E^^) 
which is smaller than the average band-gap energy {E^p) taken over the binary 
constituents: the dependence of this quantity on the ordering direction is expressed 
by the relation: < Eff < E^; 

2. Both common-anion and common-cation systems show a decreasing band-gap energy 
as the substrate lattice parameter is increased; 

3. The structures studied offer interesting opportunities for band-gap tuning as a function 
of growth conditions; the range in which the gap varies is as large as 0.7 eV in 
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(GaSb)i/(InSb)i type systems and 0.3 eV in (InAs)i/(InSb)i systems; 

4. In the common-anion (common-cation) systems the marked charge density localization 
of the CBM on the GaSb (InAs) monolayer and of the VBM on the InSb monolayer 
causes the gap to be "spatially indirect"; 

5. In the case of free standing mode elastically relaxed structures we obtain a band-gap 
value of 0.05 ± 0.05 eV in the common-anion system (semiconducting properties) and 
of -0.26 ± 0.05 eV in the common-cation system (semimetallic properties); 
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TABLES 

TABLE I. Bond-lengths (dcaSb an d di n ,sb in a.u.) and strain parameters (e GaSb and e InSb ) 
parallel and perpendicular to the growth plane for (GaSb)i/(InSb)i [111] and [001] ordered systems. 
The quantities denoted by A indicate percentage deviations from calculated bulk bond-lengths 
{dcaSb = 5.00 a.u. and dj n sb = 5.29 a.u.). 
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+3.0 % 


+1.0 % 




2.2 % 


e GaSb 


+0.029 





+0.058 


+0.029 





+0.058 


MaSb 
e _L 


-0.014 





-0.028 


-0.027 





-0.053 


dlnSb 


5.22 


5.16 


5.29 


5.25 


5.24 


5.29 


AlnSb 


-1.3 % 


-2.5 % 




-0.8 % 


-0.9 % 




e InSb 


-0.028 


-0.055 





-0.028 


-0.055 





JnSb 
e _L 


+0.016 


+0.033 





+0.030 


+0.076 
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TABLE II. Bond-lengths (di n A s and di n sb in a.u.) and strain parameters (e : and e ) 
parallel and perpendicular to the growth plane for (InAs)i/(InSb)i [111] and [001] ordered systems. 
The quantities denoted by A indicate percentage deviations from calculated bulk bond-lengths 
(dinAs = 4.96 a.u. and di n sb = 5.29 a.u.). 







[111] 






[001] 






El.Rel. 


InAs-subs. 


InSb-subs. 


El.Rel. 


InAs-subs. 


InSb-subs. 


dinAs 


5.04 


4.96 


5.12 


5.01 


4.96 


5.07 


&-InAs 


+1.6 % 




+3.2 % 


+1.0 % 




2.2 % 


JnAs 


+0.033 





+0.066 


+0.033 





+0.066 


,InAs 
e _L 


-0.019 





-0.038 


-0.036 





-0.072 


dlnSb 


5.21 


5.14 


5.29 


5.24 


5.22 


5.29 


^InSb 


-1.5 % 


-2.8 % 




-0.9 % 


-1.3 % 




JnSb 

€ \\ 


-0.031 


-0.062 





-0.031 


-0.062 





JnSb 


+0.019 


+0.037 





+0.033 


+0.078 
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TABLE III. Calculated electronic energy levels (in eV) with respect to the VBM for the [111] 
SL (neglecting s.o. coupling). The superscripts (1) and (2) indicate the two states involved in the 
repulsion mechanism. The state multiplicity is given in parentheses. 



SL State 


ZB State 


ER 


(GaSb)i/(InSb)i 
SI 


S2 


ER 


(InAs)i/(InSb)i 
S3 


S4 


r (2) 
1 lc 

r (l) 
1 lc 


Lie 

Tie 


0.50 
-0.81 


0.70 
-0.64 


-0.20 
-1.47 


0.79 
-1.02 


1.04 
-0.95 


0.36 
-1.42 


r (2) 

Tit, 

r (l) 
1 3d 


r i5D 
Ti5d 
L^v 


0(2) 
-0.11 
-1.31(2) 


0(2) 
-0.59 
-1.35(2) 


-0.42(2) 


-1.68(2) 


0(2) 
-0.21 
-1.38(2) 


0(2) 
-0.77 
-1.46(2) 


-0.18(2) 


-1.50(2) 


TABLE IV. Calculated electronic energy levels (in eV) with respect to the VBM for the [001] 
SL (neglecting s.o. coupling). The superscripts (1) and (2) indicate the two states involved in the 
repulsion mechanism. The state multiplicity is given in parentheses. 


SL State 


ZB State 


ER 


(GaSb)i/(InSb)i 
SI 


S2 


ER 


(InAs)i/(InSb)i 
S3 


S4 


r (2) 
1 lc 

p(l) 

1 lc 


Tic 


1.12 
-0.49 


1.26 
-0.43 


0.78 
-0.77 


1.46 
-0.68 


1.54 
-0.69 


1.01 
-0.96 


r (2) 

^4v 

r (l) 
1 5d 


Tl5D 
r i5D 
^5v 


0(2) 
-0.07 
-2.52(2) 


0(2) 
-0.41 
-2.52(2) 


-0.20(2) 


-2.71(2) 


0(2) 
-0.06 
-2.46(2) 


0(2) 
-0.46 
-2.46(2) 


-0.25(2) 


-2.70(2) 
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TABLE V. Band-gap energies (in eV) for (GaSb)i/(InSb)i [111] systems and (InAs)i/(InSb)i 
[111] obtained from unperturbed LDA calculations (Eg2p), with the introduction of spin-orbit 
coupling {EgJ^ A ) and corrected starting from experimental data (£*™ p ). We also report the 
calculated bowing parameter (bj^ ) for the different systems considered. 





ER 


(GaSb)i/(InSb)i 
SI 


S2 


ER 


(InAs)i/(InSb)i 
S3 


S4 


jpunp 

gap 


-0.81 


-0.64 


-1.47 


-1.02 


-0.95 


-1.42 


rpLDA 
gap 


-1.05 


-0.90 


-1.59 


-1.26 


-1.20 


-1.49 


rpemp 
gap 


0.05 


0.20 


-0.47 


-0.26 


-0.21 


-0.51 


6 [1H] 
u gap 


1.92 


1.32 


4.08 


2.42 


2.18 


3.34 



TABLE VI. Band-gap energies (in eV) for (GaSb)i/(InSb)i [001] systems and (InAs)i/(InSb)i 
[001] obtained from unperturbed LDA calculation (E^p), with the introduction of spin-orbit 
coupling {Eg^ A ) and corrected starting from experimental data (Eg™ p p ). We also report the 
calculated bowing parameter for the different systems considered. 





(GaSb)i/(InSb)i 
ER SI S2 


(InAs)i/(InSb)i 
ER S3 S4 


Tpunp 
gap 

rpLDA 
gap 

jpemp 
gap 


-0.49 -0.43 -0.77 
-0.74 -0.69 -0.96 
0.33 0.41 0.13 


-0.68 -0.69 -0.96 
-0.88 -0.89 -1.07 
0.11 0.10 -0.09 


J001] 
"gap 


0.68 0.48 1.56 


0.90 0.94 1.66 
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TABLE VII. Band-gap energies E^ A , spin-orbit coupling included, and crystal- field splittings 
Acf (in eV) for common-anion systems. 





(GaSb)i/(InSb)i 
ER SI S2 


(GaSb)i/(GaSb)i 
ER SI S2 


(InSb)i/(InSb)i 
ER SI S2 


rpLDA 
gap 

A CF 


-1.05 -0.90 -1.59 
0.11 0.59 -0.42 


-1.15 -1.14 -1.66 
0.19 0.58 -0.42 


-0.34 -0.26 -0.88 
0.14 0.57 -0.50 


TABLE VIII. Band-gap energies Eg^ A , spin-orbit coupling included, and crystal-field 
splittings Acf ( m eV) for common-cation systems. 




(InAs)i/(InSb)i 
ER S3 S4 


(InAs)i/(InAs)i 
ER S3 S4 


(InSb)i/(InSb)i 
ER S3 S4 


rpLDA 
gap 

A CF 


-1.26 -1.20 -1.49 
0.21 0.77 -0.18 


-1.26 -1.24 -1.44 
0.14 0.57 -0.39 


-0.29 -0.25 -1.09 
0.11 0.68 -0.47 
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TABLE IX. Angular decomposition relative to the muffin tin charge density (for s (Q s ) and p 
(Qp) components) of the different states, neglecting s.o. coupling, for the (GaSb)i/(InSb)i [111] 
systems. 



State 


Ga 


SI 
Sb 


In 


Sb 


Ga 


ER 

Sb 


In 


Sb 


Ga 


S2 
Sb 


In 


Sb 




Qs 
Qp 


0.065 
0.037 


0.057 
0.187 


0.006 
0.060 


0.005 
0.134 


0.003 
0.040 


0.001 
0.218 


0.014 
0.048 


0.017 
0.184 


0.007 
0.074 


0.001 
0.115 


0.045 
0.014 


0.023 
0.241 


r (2) 


Q s 
Qp 


0.000 
0.023 


0.000 
0.122 


0.000 
0.061 


0.000 
0.321 


0.000 
0.025 


0.000 
0.117 


0.000 
0.063 


0.000 
0.306 


0.000 
0.026 


0.000 
0.114 


0.000 
0.063 


0.000 
0.294 


r (l) 
1 lc 


Qs 
Qp 


0.262 
0.019 


0.177 
0.027 


0.074 
0.002 


0.072 
0.063 


0.345 
0.005 


0.232 
0.019 


0.047 
0.006 


0.045 
0.001 


0.349 
0.011 


0.217 
0.031 


0.031 
0.005 


0.030 
0.001 


r (2) 
1 lc 


Qs 
Qp 


0.028 
0.073 


0.008 
0.024 


0.188 
0.000 


0.161 
0.025 


0.008 
0.057 


0.000 
0.006 


0.228 
0.001 


0.180 
0.052 


0.000 
0.001 


0.003 
0.106 


0.223 
0.033 


0.180 
0.004 
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FIGURES 

FIG. 1. Calculated highest valence band energy levels (Ei, E2 and E3) and lowest conduction 
state (E c ) at T versus substrate lattice constant for the [111] (a) common-anion and (b) 
common-cation SLs considered. The center of gravity of the topmost valence bands is taken as 
zero. 

FIG. 2. LDA band-gap (E LDA - solid line) and corrected band-gap (E emp - dashed line) as a 
function of the substrate lattice constant for the [111] and [001] (a) (GaSb)i/(InSb)i SLs and (b) 
(InAs)i/(InSb)i SLs. Filled (empty) circles indicate the average value (E ave ) of the experimental 
(LDA) band-gap in the pure binaries. 

FIG. 3. Charge density distribution (in units of 0.5 e/unit cell) for the (a) T\ v state and for 

(2) 

the (b) T 3v state in the elastically relaxed [111] ordered common-anion structure. 

FIG. 4. Charge density distribution (in units of 0.5 e/unit cell) for the (a) T^} (CBM) state 
and for the (b) T^} state in the elastically relaxed [111] ordered common-anion structure. 
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